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The interaction between the radiation field in a microcavity with a mirror undergoing damping 
oscillation is investigated. Under the heavily damping cases, the mirror variables are adiabatically 
eliminated. The the stationary conditions of the system are discussed. The small fluctuation ap- 
proximation around steady values is applied to analysis the antibunching effect of the cavity field. 
The antibunching condition is given under two limit cases. 
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' 1. Introduction 

o ■ 

There has been interest in the investigation of the non-classical behavior of the light field, for example, the squeezed 
] state of the electromagnetic field [0-|4|. A single mode or a multimode electromagnetic field is squeezed when the 
fluctuation of its one quadrature component is reduced to below the standard quantum limit. But the fluctuation of 
^ i' its conjugate partner will be enlarged in terms of the uncertainty relation. The squeezed electromagnetic field could 

be applied to the optical communication and ultrasensitive gravitational wave detection. 
^ [ Recently, it has been shown that a cavity with free oscillating mirror might be employed as a model for squeezing. It 
was firstly proposed by Stenholm without taking into account the effect of the fluctuation on the oscillating mirror . 
Since then, this model is generalized to the coupling of the system to the external world ■ But up to now, another 
non-classical behavior of the radiation field in a microcavity with a movable mirror , such as antibunching effect, isn't 
still studied. 

In this paper, we will discuss the antibunching effect of the radiation field in a microcavity with a mirror undergoing 
^ ' heavily damping oscillation. In section 2, we firstly give a quantizing Hamiltonian including external dissipative effect. 
The master equation of the reduced density matrix for the system is given. We convert the master equation into a c- 
number equation (the Fokker-Planck equation) by using the two-mode positive P representation. Then the stochastic 
equation corresponding to the Fokker-Planck equation is obtained. In section 3, we discuss the stability of the system 
under the good cavity limit and linearize the stochastic equation around stable values. In section 4, antibunching 
effect is investigated and antibunching condition of the system is given. Finally we give the conclusion of this paper. 



2. Hamiltonian and Master equation 

We firstly consider the interaction between a single mode cavity field and a movable mirror. Based on reference 
we have a effective quantizing Hamiltonian 

Hq = hujcd^a + huj,nb+b - nga+a{b+ + b) + ih{E{t)a+ - E*{t)a) (1) 

where E{t) is proportional to the amplitude of the external driving field of the cavity mode. E*{t) is complex conjugate 
of E{t). ujc is the cavity frequency. a{a^) are annihilate (create) operators of the cavity field, ujm is a frequency of 
the mirror. &(&+) are annihilate (create) operators of the mirror, g — ( 2rnuj ) • '^^ is a mass of the mirror and L is 
the equilibrium cavity length. 



If the cavity is bad, and the mirror is damped by the circumstance when it is moving, then the Hamiltonian (1) 
need correct to add the dassipative effect of the circumstance. So we have 

H = Ho + ha+Ti + har+ + hb+r2 + hbr+ (2) 

where ri(r5^) is the reservoir operators of the cavity. We also model the damping effect of the mirror as the result of 
the interaction between the mirror and the many harmonic oscillator. They satisfy the Markovian correlation function 
(where for simplicity, we only consider the zero temperature case): 

<T,{t)T+{t') > =j,5{t-t') 

<T+(t)r,{t')> =0 (3) 

with i = 1, 2. 7i is the decay rate of the cavity field. 72 is the damping coefficient of the motion mirror. The master 
equation of the reduced density matrix for the system is 

dp 1 

— = — [Ho,p] + 7i(2apa+ - pa+a - a+ap) 

+ -f2{2bpb+ - pb+b-b+bp) (4) 

We could convert the operator eq. (4) into a c- number equation (the Fokker-Planck Equation) by using the two- mode 
positive P representation . This representation ensures that the P function exists as a well-behaved distribution 
which is singular when the usual Glauber-Sudarshan P representation was applied. That is: 

dP d 

— = [- — (iWcOi - igaio^ - igaia2 + 71"! - E{t)) 
at oai 

d 



- — {iuj,na2 - igajai +72^2) 
oa2 

d_ 

Z+ 



where ai and a^, and a2 and a J are no longer complex conjugate of each other, instead they are independent 
complex variables. In terms of the Ito ruler, the Fokker-Planck eq.(5) is equivalent to the following set of the stochastic 
equations. 



9a 1 

da+ 



-iujcCti + igaia2 + iga\a2 — "fiO-i -f E{t) + (6a) 
'^cCf-i — iga^a^ — iga^a2 — ^lO.^ + E*{t) + (6b) 



= -iu!ma2 + igajai - 7202 + (6c) 



iuijna^ - iga^ai - 720;^ + (6d) 



dt 

where Fq.^ and F^+ are Gaussian random variables with zero mean. Their correlation functions satisfy: 

< r„, (i)r„, it') > = < r,, (t)r„, (f) iga^Sit - t') (7a) 

< r„+ (t)r„+ [f] > = < r„j (t)r„+ (f) ~iga+6it - 1') (7b) 

< F„, (t)r^+ (f) > - < (<)r„, (f) >^ o (7c) 

with Z = 1 or 2 and m 1 or 2. The above correlation functions could be obtained from the matrix elements of the 
diffusion matrix of eq.(5) by using Ito ruler. 

3. Adiabatic elimination of the mirror variable and stability analysis 
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In the good-cavity limit, that is, the damping cocfRcicnt of the mirror is much larger than the decay rate of the cavity 
field (72 J> 71). "This means that the decay time of the cavity field is much longer than the decay time of the 
mirror. So the mirror variables could be adiabatically eliminated from eqs.(6a-6b) . Firstly we find a steady-state of 
the mirror variables . In the case of the steady state, 0:2 = and a2 = 0. So we have from eqs.(6c-6d) 



Prom eqs.(8a-8b), we obtain: 



- iuima2 + igc(t<^i - 72Q!2 + T^^ = (8a) 
iuimCtt - igatai - 720;^ + r^+ = (8b) 



a, = ^"-+^^^^"^ (9a) 

72 + lUJrn 

^ r + - igaia^ 
at = : (9b) 

72 - lUJm 



We substitute eqs.(9a-9b) into eqs.(6a-6b) and eliminate the mirror variables. Then we have: 



= -iuicai + iG{u>rn)aiai - licti + T (10a) 



dt 



= itocaj - iG{ujra)aJ^ai - 7ia]^ + r+ (10b) 



with G(a;„) = and 



igax tgax 

1 = — —. — H : — r + -i-rai (iia) 

72 + 72 - lU)m ^ 

r+ = -^?^r., - ^^r„. + (iib) 

72 + «Wm 72 - "2 "1 



Using eqs.(7a-7c) and eqs.(lla-llb), we calculate the correlation function 



< r(i)r+(i') > = < r+(t)r+(t) >= ^J'^'" 8{t - 1') (i2a) 

<T+(t)T+(t')> = ^—^6(t-t') (12b) 

72 - «Wm 

< r(t)r(t') > = — 5{t - 1') (12c) 

72 + lUJm 

The eqs.(lOa-lOb) are difiicult to be solved. In general, we are interested in the properties of the steady state. So 
we denote the steady values of ai and by ao and respectively. We assume that the system has a small 
approximation around the steady values, namely 

ai{t) = ao + 5ai{t) , . 

at{t) = at + 6at{t) ^'^^ 

So non-linear eqs.(lOa-lOb) are simplified into the following linear equations around the steady values. 

— ( '^"1 ^ ^ ( "^"^c - 71 + i2G{(jJm)C(oC(0 iG{u)m)cto \ ( 

dt V ^oii J ~ \ -iG{ijJm)oit'^ ioJc - 71 + i2G{ijJm)atao J \ ^ol^ 



-:?rx — 



+ I _ Vo+l ( n^i^t) ' ^^^^ 

where n = ajao and r\i(t) statisfy delta correlation function: 

<r]i{t)> =0 (15a) 

<Vi{t)Vj{t') > =6^j6{t-t') (15b) 
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We abbreviate eq.(8) as following: 

^^Sa{t)^^Aa{t) + Dirf{t) (16) 

But one of the important feature of eq.(14) is whether the steady solutions are stable, that is, when ai and 
somewhat deviate from ao and , whether they will still return to steady values. So we neglect the fluctuation forces 
of eq.(16) and have: 

pd{t) = -Ad{t) (17) 

In order to investigate the stationary of the system, we seek the solutions of the form e~^* of eq.(17). The eigenvalues 
A are determined by the equation 

\A-XI\=0 (18) 
where / is a identity matrix. The solutions of the eq.(18) are 

A = 71 ± \/ {ujc - 3G{LLijn)n){uJc - G{oJm)n) (19) 

This equation indicates if G{u!m)n < loc < 3G{L0m)n-, the real part of the eigen-solution of the eq.(18) is positive. 
The system is stable. When lUc < nGluJm) or lOc > 3nG{uJm), only if 71 > ^y {ujc — iG{uJm)n){uJc — G{ijJm)n) then the 
system is also stable. So only we chose the proper parameter, the system always may reach to stability. The small 
fluctuation approximation is appreciate. 

4. Second-order correlation function 

The intensity correlation is another quantity of the experimental interest besides the first order optical coherence. 
It's truly photon correlation measurement. Theoretically, people have defined a second-order correlation function 
to investigate the joint photocount probability of detecting the arrival of a photon at one time and other ptoton at 
another time. For zero time delay. The second order correlation function is ||ll[] 

g^-^(o)^ <;Yy (20) 

Under the positive P representation we keep the second order terms of then eq(20) becomes into: 

5' ' == T r = l + 5 + T? +4 21 

< a+a >2 a+a 

In order to calculate the correlation functions < {6a)'^ > et.al, we need calculate matrix: 

Q _ f < 5(\5a > < SaSa^ > \ .^^ 
~ \ < 5a+5a > < 6a+Sa+ > J ^ ' 

The matrix C could be obtained by calculating of Sai and 5ai from the linearized eq.(14). It is also could be given 
by 0: 

^ D.DetjA) + {A- ITr{A))D{A ITr{A)Y 

2Tr{A)Det{A) ' ^ ' 



after tedious calculation, we have: 



al[2M-fl + 2NG[B + ^71) - i2jiBM - G^n'^iM + M* 
471(72 + 52 _G2n2) 



Gil - O22 - ... I R2 7^7Zt{ i^4aj 



„ _„ _ 2jV(72 + B^) + ^7iGn2(A/* - Af ) - GBn^jM + M*) 

- ' 471(7?+ 52 -G2n2) ^^^^^ 



with 
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M ^ 2£ — N = B = ujc- 2Gn (25) 

So the second order correlation function is: 

_^ ^ G[{ujc - 2nG){2nj2(^c - ^^726' - 7it^„0 + n(7?72 + ^^^^^72 - 27itj,„G'n)] 
^ njiujmhf - {c^c + 5nG){uJc - nG)] 

From this equation, we see that the antibunching of the cavity field appear when the second term of the eq.(26) is 
a negative number. This condition could be satisfied after we chose some proper parameters. Now we consider two 
limit cases. In the case of the strong cavity field, we only keep the terms including n^. The second term of eq.(26) 
is positive, no antibunching appears. But under case of the weak cavity field, namely when n ^ 0, the cavity field 
presents antibunching behavior. 

5. Conclusion 

The interaction between the radiation field in a microcavity with a movable mirror undergoing heavily damping 
oscillation is investigated. Under the heavily damping cases, the mirror variables are adiabatically eliminated. The 
stationary condition of the system is given. The small fluctuation approximation around steady values is applied to 
analysis the antibunching behavior of the cavity field. In the case of the strong cavity field, no antibunching appears. 
But under the case of the weak cavity field, cavity field presents the antibunching effect. 
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